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Abstract 

Within the "lowest Landau level approximation", we develop a method to 
find the ground state of a 2d system of interacting particles confined by a 
parabolic potential. 
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The problem of finding exact eigenstates of a 2d system of particles interacting by the 
two-body central forces V{r) and confined by an external parabolic potential arises in many 
subfields, including cooled and trapped fermionic atoms, Bose-Einstein condensates in traps, 
quantum dots and electrons in magnetic field. Assuming spherical symmetry, we are inter- 
ested here in the lowest energy states of such system as a function of exact quantum numbers. 
The lowest Landau level approximation is equivalent to the weak interaction limit 

V{r) < huj 

where u is the oscillator frequency of the external trapping potential. The effective Hamilto- 
nian of the A^-body system is the sum of the energy of degenerate level with total angular 
momentum L (/i=m=l) 

" N 



H = uj{L + N)+V = P{L) 



P{L) (1) 



which comes from the noninteracting Hamiltonian in the 2d parabolic trap, and the projected 
interaction V. Here, P{L) is the projector onto the space of lowest Landau level, S. The 
space 5* is formed by the homogeneous polynomials in Zi=Xi+iyi of degree L times the 
Gaussian factor |0)=e^^ ^ Hereafter, we set uj=l for simplicity. The polynomials are 
antisymmetric (symmetric) for identical fermions (bosons). The eigenstates of the H (|I]) are 
found from diagonalization of V in the space S. 

The effective interaction V can be found in the form p| 

lL/2] 

l> = E V,^ S2kiB'' - B''^-') + (.2fc+2 - S2k+,)B''^\ (2) 

N 

where B'' = I] (a| — a])'^(aj — aj)'' and a] = Zi/2 — d/dz* and Oj = z*/2 + d/dzi are the 

j>i 

standard ladder operators. Here, B^=N{N~l)/2 and B^^=0; [k] denotes the integer part 
of k. In (I), Sk=J^§M{k+l,l-t) V{^/2i) with M the Kummer function. 

In order to find the ground state without solving the whole spectrum, we use the following 
idea 0. Suppose that the Hamiltonian can be written as a sum 

V = Vo + Vs, (l) V^5|0) = 0, (it) Vo\0) = ErninlO), (3) 



such that the second term (i) annihilates a state |0), and (ii) |0) is the ground state for the 
first term (this ground state may be degenerate). If (iii) the second term is non-negative 
definite (i.e., it has no negative eigenvalues) 

Vs > 0, (4) 

then the state |0) will still be the ground state of the full Hamiltonian Vq + Vs, with the 
same eigenvalue Emin HI- 

In the case of Bose system, the polynomials of the basis states are Sl = PsttllPs, where 
TTiL = J'y^2---^n is the monomial with a given partition of the integer L = J^l-n, and Ps 

n 

is the operator of symmetrization over variables. The above decomposition (0,§) with 
the properties (i) and (ii) can be found by inspecting action of (0) on the basis states. The 
simplest basis state 

\0l) = PsZiZ2...zl\0), (V - Vo)\Ol) = 0, (5) 
is annihilated by any term in (|^), except for Vq. The term Vq can be written as 

Vo = {N/2)[{N - l)so - (L - v){s, - S2)]. (6) 

Here, 

1 ^ 

is the collective contribution to the angular momentum, which is also a conserved quantum 
number, [v, V] = 0. The allowed values of v are 0, 1, 2, .., L — 2, L [0]. The state (^) can be 
represented as a sum 

N L-v 

\0l) = E \^L,v), \Ol,v) = Z-PsUi^k - Z)\0), (7) 

v=0 k=l 

where Z=J2k=i^k/ N is the collective variable. Each term 10^,^,) is the eigenvector of v 
and therefore of Vq. Each f-sector is represented by a single state 10^^^,) in expansion (|^. 
The spectrum of Vq consists of L equidistant degenerate levels. Conditions (i) and (ii) are 



therefore satisfied for the states \Ol,v)- In particular, this means that those are eigenstates 
of V with eigenvalues given by (^) 0, see also [Q. 

If the term Vs = V — Vq is non-negative definite, the condition (iii) is satisfied and each 
state \Ol,v) will be the ground state in corresponding sector L,v. In fact, it is possible to 
find a broad universality class of the potentials V{r) which give Vs = V — Vq > 0. We work 
in the space S and use permutation symmetry to write PsVsPs in the form 

Ps E VS,^,PS = ^^^^ ^^ PsVs,12Ps = ^^^~^^ PsPl2Vs,uPl2Ps (8) 

i>j 

where Vs,ij is the contribution from pair of particles i,j to Vs (0,3), and pij is the operator 
of symmetrization over this pair. The condition puVs^uPu > is a sufficient condition for 
PsVsPs > 0. The eigenvalues of Pi2Ks',i2Pi2 can be easily evaluated, using the full space 
of monomials m^,. As a result, the condition pi2Vs^i2Pi2 > reduces to the following set of 
integral inequalities 



j dtV{ 



r)e-^'/' 



An / j,4' 



22" (2n) 



-l+n\l 



> 0, (9) 





which must be satisfied for integer n, < 2n < L. These conditions describe the broad 
universality class of interaction potentials V{r), for which the lowest energy states in sec- 
tors with quantum numbers v and L<N are given by {Ol^^) [see Eq.(^], and their ener- 
gies are equal to the eigenvalues of Vq, Eq.(^, with so=J^dte~*V{\/2i), Si—S2=f^dte~*' 
— ^^V{\/2i). Conditions ® are valid, in particular, for the Gaussian interaction V{r) = 

— / B?' _^ 

Uq '^ ^j^2 ,Uo>0, for the (5-function interaction V=Uo S{r),Uo>0 and for the 2d log-Coulomb 
interaction V=Uolog{l/r), Uq>0 and for many others potentials, which are predominantly 
repulsive. These results are obtained for the Bose case. Generalization of the method for 
the Fermi case is possible. 
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